Gravitational lensing in fourth order gravity 
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Gravitational lensing is investigated in the weak field limit of fourth order gravity in which the 
Lagrangian of the gravitational field is modified by replacing the Ricci scalar curvature R with an 
analytical expression f{R). Considering the case of a pointlike lens, we study the behaviour of the 
defiection angle in the case of power -law Lagrangians, i.e. with f{R) oc R". In order to investigate 
possible detectable signatures, the position of the Einstein ring and the solutions of the lens equation 
are evaluated considering the change with respect to the standard case. Effects on the amplification 
of the images and the Paczynski curve in microlensing experiments are also estimated. 
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I. INTRODUCTION 



The last decade observational results, including 
Type la Supernovae WLCOsmic microwave background 
anisotropy spectrum [2 13 and large scale structure jj| , 
have radically changed our view of the universe leading 
to revise the old cosmological standard model. On the 
contrary, the nowadays standard scenario, referred to as 
the concordance model, although still spatially flat, as- 
sumes that the visible baryonic matter constitutes only 
~ 4% of the matter - energy content, while the most of 
the budget is in the form of dark components. This 
dark sector is assumed to be constituted for one third by 
dark matter, characterized by ordinary thermodynamical 
properties and mainly clustered in galaxies and clusters 
of galaxies, and for two third by a negative pressure fluid 
dubbed dark energy, which drives the large scale cosmic 
acceleration and does not cluster at large and small scales 
hence not influencing galactic dynamics. 

Investigating the nature and the properties of the dark 
side of the universe represents the most daunting and yet 
fascinating challenge of modern cosmology so that it is 
not surprising that there are a lot of theoretical propos- 
als on the ground. Nevertheless, both in the case of the 
old problem of dark matter and the more recent issue of 
dark energy, few approaches seem to furnish a quite nat- 
ural and coherent framework within which to explain the 
observational results. Still more rare are models which 
provide a unified view of both dark matter and dark en- 
ergy as a two separate manifestations of a single scale - 
dependent gravitational phenomenon. 

From this point of view, higher order theories of gra vity 
(both in the metric |1H|1||| and the Palatini [9lll0l[Tl| 

formulations) represent an interesting approach able to 
fruitfully cope with both dark matter and dark energy 



problems^. On one hand, it has been widely demon- 
strated that such theories can agree with the cosmologi- 
cal observations of the Hubble flow jlj, [3 and the large 
scale structure evolution [15i] . On the other hand, in the 
weak field limit , the g ravitational potential turns out to 
be modified |lall7llla . ll9ll20 | in such a way that interest- 
ing consequences on galactic dynamics may be achieved 
without violating, at the same time, the constraints on 
the PPN parameters ^3 from Solar System tests. 

It is worth remembering that one of the first experi- 
mental confirmations of Einteinian general relativity was 
the deflection of light observed during the Solar eclipse 
of 1919. Since then, gravitational lensing (i.e., the de- 
flection of light rays crossing the gravitational field of a 
compact object referred to as lens) has become one of 
the most astonishing successes of general relativity and 
it represents nowadays a powerful tool able to put con- 
straints on different scales, from stars to galaxies and 
cluster of galaxies, to the large scale structure of the uni- 
verse and cosmological parameters (see ^^ and [23 for 
comprehensive textbooks). 

Being intimately related to the underlying theory of 
gravity in its Einstein formulation, it is quite obvious that 
modifying the Lagrangian of the gravitational field also 
affects the theory of gravitational lensing. It is therefore 
mandatory to investigate how gravitational lensing works 
in the framework of higher order theories of gravity. On 
the one hand, one has to verify that the phenomenol- 
ogy of gravitational lensing is preserved in order to not 
contradict those observational results that do agree with 
the predictions of the standard^ theory of lensing. On 
the other hand, it is worth exploring whether deviations 
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An alternative yet similar approach is the one based on scalar - 
tensor theories of gravity (see, e.g., 1121 and references therein). 
In the remaining of the paper, we will refer to the theory of grav- 
itational lensing based on Einstein general relativity as the stan- 
dard or Einsteinian theory, while as corrected theory we mean 
the one based on fourth order gravity. 



from classical results of the main lensing quantities could 
be somewhat detected and work as clear signatures of a 
modified theory of gravity. As a first step towards such 
an ambitious task, we consider here power -law fourth 
order theories, i.e. we replace the Ricci scalar R in the 
gravity Lagrangian with the function f{R) oc i?", and 
investigate how this affects the gravitational lensing in 
the case of a pointlike lens. 

The paper is organized as follows. After a short sum- 
mary of f{R) theories, we derive, in Sect. Ill, the general 
expression for the defiection angle of a pointlike lens in 
the weak field regime. In particular, we discuss why this 
result holds whatever is the theory of gravity which en- 
ters in determining the expression for the gravitational 
field. Sect. IV is devoted to evaluating the defiection an- 
gle and deriving the lens equation for the case of the f{R) 
theory we are considering. Since our results refer to the 
case of a pointlike lens, we may compare and contrast 
them with the classical ones in the microlensing regime. 
This issue is extensively discussed in Sect.V where we 
study the deviations from the classical results regarding 
the position of the images, the Einstein angle, the to- 
tal amplification and the Paczynski lightcurve, while, in 
Sect. VI, we summarize and conclude. 



II. BASICS OF f{R) THEORIES 

Higher order theories of gravity, also referred to as 
f{R) theories, represent the most natural generalization 
of the Einsteinian general relativity. To this aim, one 
considers the gravity action : 



Eq.(|2l and adopting the Robertson - Walker metric, it is 
possible to show that the Friedmann equations may still 
be written in the usual form provided that 'ineffective 
curvature fluid is added to the matter term with energy 
density pcurv and pressure p^urv = WcurvPcurv depending 
on the choice of f{R) as : 



f'{R) \2 



[f{R)-Rf'{R)]~iHRf"{R) 



= -1 



Rf"{R) -f R Rf"'{R) - Hf"{R) 



[f{R) - Rf'{R)] /2 - 3HRf-{R) 



(3) 



(4) 



As a particular case, we consider power -law f{R) theo- 
ries, i.e. we set : 



/(i?) = /oi?" 



(5) 



with n the slope of the gravity Lagrangian (n = 1 being 
the Einstein theory) and /g a constant with the dimen- 
sions chosen in such a way to give f{R) the right physical 
dimensions. It has been shown that such a choice leads 
to matter only models able to fit well the Hubble diagram 
of Type la Supernovae without the need of dark energy 
[li^ and could also be reconciled with the constraints on 
the PPN parameters |21|. 



III. THE DEFLECTION ANGLE 



A 



[fiR)+Cm] 



(1) 



where f{R) is a generic function of the Ricci scalar cur- 
vature R differentiable at least up to the second order 
and Cm is the standard matter Lagrangian. The choice 
f{R) = R + 2A gives the general relativity including the 
contribution of the cosmological constant A. Varying the 
action with respect to the metric components (/^^, one 
gets the generalized Einstein equations that can be more 
expressively recast as p : 



Although generalizing the gravity Lagrangian has a 
deep impact on both the cosmology and the local dy- 
namics, it is nevertheless easy to understand that the 
derivation of the lens equation is formally the same as 
that in standard general relativity. Indeed, the basic as- 
sumption in deriving the lens equation is that the grav- 
itational field is weak and stationary. In this case, the 
spacetime metric reads : 



ds' 



2$ 



'dt' - 1 



2$ 

72" 



Sijdx'dx^ (6) 



G,. = -^UffM^[/(i?)-i?/'(i?)]+/'(i?);M. 



9^^0f{R) 



IW) 



(2) 



where G^v — Rf^i, — {R/2)g^i, and T)i" are the Einstein 
tensor and the standard matter stress - energy tensor re- 
spectively and the prime denotes derivative with respect 
to R. The two terms f'{R).^^ and Uf'{R) imply fourth 
order derivatives of the metric g^n, so that these models 
are also referred to as fourth order gravity. Starting from 



where $ is the gravitational potential and, as usual, we 
have neglected the gravitomagnetic term'^. Since light 
rays move along the geodetics of the metric ©, the lens 
equation may be simply derived by solving ds^ = and 
keeping only terms up to second order in v /c. It is worth 
stressing that such a derivation holds whatever is the 



^ Gravitomagnetic effects originate from the mass current in the 
lens and, under some circumstances, could give rise to effects 
that could be detected with future interferometric satellites (see, 
6-g-i l24l and references therein). 



theory of gravity provided that one can still write Eq. © 
in the approximation of weak and stationary fields. As 
a fundamental consequence, the lensing deflection angle 
will be given by the same formal expression that holds in 
general relativity ^22, 23i| 



Vi $d/ 



(7) 



with : 



V 



V - e(e-V) 



(8) 



where e is the spatial vector^ tangent to the direction 
of the light ray and dl — y/6ijdx^dx^ is the Euclidean 
line element. The integral in Eq.Q should be performed 
along the light ray trajectory which is a priori unknown. 
However, for weak gravitational fields and small deflec- 
tion angles, one may integrate along the unperturbed di- 
rection. In this case, we may set the position along the 
light ray as : 



^ + le 



(9) 



with ^ orthogonal to the light ray. Assuming the poten- 
tial only depending on r = |r| = (^^ + Py/^ (as for a 
pointlike or a spherically symmetric lens), the deflection 
angle reduces to : 



a 



c^ ./„^ \r dr J 



n2 



(10) 



where we have assumed that the geometric optics approx- 
imation holds, the light rays are paraxial and propagate 
from infinite distance. 

Ea. HlU|) is general and holds whatever is the under- 
lying theory of gravity. This nice result does not mean 
that the defiection angle is the same as in the Einsteinian 
general relativity. Indeed, the link with the underlying 
gravity theory is represented by the gravitational poten- 
tial $ which is determined by solving the corresponding 
generalized Einstein equations for the weak field metric 
©. As a result, <I> can be no more the standard Newto- 
nian potential and deviations from the Keplerian scaling 
1/r enter the game leading to interesting consequences. 
From the lensing point of view, this means that, although 
Ea. (|10|l still holds, the deflection angle differs from the 
one of general relativity in a way that critically depends 
on the expression chosen for the function /(i?) entering 
the gravity Lagrangian. 



* Here, quantities with an over - arrow are vectors, while the versor 
will be denoted by an over -hat. 



Before deriving explicitly the deflection angle in f{R) 
theories, it is worth discussing an important point. The 
above considerations rely on the assumption that the 
weak field metric may be still be written as in Eq.@ 
in the case of alternative theories of gravity. That is, 
one is implicitly assuming that the solution of the field 
equations in the low energy limit shares the same for- 
mal structure as the Schwarzschild solution which Eq. @ 
reduces to when ^(r) is the Newtonian 1/r potential. Ac- 
tually, under the hypothesis of weak gravitational fields 
and slow motions (both well verified in the astrophysical 
phenomena of interest), by virtue of the Birkhoff- Jensen 
theorem, the weak field metric is : 



ds^ = A{r)dt^ - B{r)dr^ - r'^dVt^ 



(11) 



de^ 



is the line element on the 



where dVi^ 

unit sphere. Inserting Ec ljllll into the Einstein field 
equations one obtains the classical Schwarzschild solu- 
tion by setting A{r) = 1/ B{r) = 1 + 2^{r)/c^ and using 
the limit <3?/c^ << 1 with <I>(r) the Newtonian potential. 
One could argue that this is no more the case when al- 
ternative theories of gravity are considered so that one 
should take A{r) ^ 1/ B{r) and then deriving the lens 
equation in the general PPN formalism as in |25j. How- 
ever, considering higher order theories of gravity leads to 
higher order field equations and hence widens the set of 
solutions. It is therefore possible, a priori, to find out 
solutions of the low energy limit of the field equations 
under the assumption A{r) = 1/B{r) so that the space- 
time metric indeed reads as in Eq.® and the deflection 
angle is given by Ea. (|10|l also for f{R) theories. 



IV. THE POINTLIKE LENS 

Eq. H1U|I allows to evaluate the deflection angle provided 
that the source mass distribution and the theory of grav- 
ity have been assigned so that one may determine the 
gravitational potential. As a first step, we consider here 
the case of the pointlike lens. Note that, although be- 
ing the simplest one, the pointlike lens is the standard 
approximation for stellar lenses in microlensing applica- 
tions 22, 2q. Moreover, since in the weak field limit a 
is an additive quantity, the defiection angle for an ex- 
tended lens may be computed integrating the pointlike 
result weighted by the deflector surface mass distribu- 
tion under the approximation of thin lens 22] (i.e., the 
mass distribution extends over a scale which is far smaller 
than the distances between observer, lens and source). 

Given the symmetry of the problem, it is clear that 
we may deal with the magnitude of the deflection angle 
and of the other quantities of interest rather than with 
vectors. In the approximation of small deflection angles, 
simple geometrical considerations allow to write the lens 
equation as : 






(12) 



which gives the position 6 in the lens plane of the images 
of the source situated at the position 9s on the source 
plane^. Note that the lens and the source planes are de- 
fined as the planes orthogonal to the optical axis which 
is the line joining the observer and the centre of the lens. 
Here and in the following, Z?;, Z?s, Dis are the angu- 
lar diameter distances between observer -lens, observer - 
source and lens - source respectively. 

In order to evaluate the deflection angle, we need an 
explicit expression for the gravitational potential $ gen- 
erated by a pointlike mass ni. This has been derived in 
detail in |23 for the case of power -law f{R) theories. 
Solving the low energy limit of the field equations (HJ in 

vacuum (T^™ — 0) for f{R) given by Eq.lO), one obtains 
for the gravitational potential : 



^{■r)=-—- 
2r 



l+(- 



(13) 



with: 



/3 



12n2 



7n - 1 - V36n4 _^ i2n^ - 83n2 + 50n + 1 



6n2 - 4n + 2 



(14) 

For n = 1, (3 — and the potential reduces to the Newto- 
nian one as expected. While the slope P of the correction 
term is a universal quantity (since it depends on the expo- 
nent n entering the gravity Lagrangian), the scalelength 
Tj, is related to one of the integration constant that have 
to be set to solve the fourth order differential equations 
of the theory. As such, re is expected to be related to 
the peculiarities of each system and can therefore take 
different values depending on its mass and typical scale. 
Inserting Ea. (|13|l into H1U|) and using ^ = Di6, after 
some algebra, we finally get for the deflection angle : 



a 



2Gm f ^ 



1 



^(1 _ p)Til - (i/2) ( e 



2r(3/2-/3/2) 

(15) 
which is defined^ only for < /? < 2. In the follow- 
ing, however, we will only consider the narrowest range 
< /3 < 1 since, as discussed in |20|, in this case, the 
modified potential offers the possibility to fit galactic ro- 
tation curves without dark matter. As a cross check, note 
that, for n = 1, it is /3 = and : 



lim a(^, /?) = UGR 



AGm 



so that the classical^ result of Einsteinian relativity is 
recovered. On the contrary, for /3 = 1, we get a = aoR/'i 
so that the net effect of the correction term is to decrease 
the deflection angle with respect to the classical one. 

As we have quoted above and explicitly shown in [20|, 
the modified gravitational potential leads to an increase 
of galactic rotation curve acting as a sort of effective dark 
matter filling the gap between the Newtonian and the ob- 
served rotation curve. Given that the classical deflection 
angle is proportional to the mass of the system, one could 
expect that the modified potential leads to an increase 
of the deflection angle while we get the opposite result. 
However, there is an important difference explaining this 
counterintuitive behaviour. Indeed, the deflection angle 
depends on the integral of X = (1 / r)d<^ / dr . Comparing 
the Newtonian potential with the modified one, we get : 



$(r) 
•^GRir) 



lir) 1 



iGRir) 2 



l + (l-/3) 



so that 



^r) > <^GRir) ^=> r>rc 



Ar) > lGR{r) ^^ r > (1 - pyf'r, . 

In order an increase of the rotation curve (and hence a 
good fit without dark matter), the first condition should 
be met which is easy to realize on galactic scales given 
typical values of Tc- On the contrary, for a pointmass lens, 
since /3 < 1, the integral entering the deflection angle 
T{r) is actually always smaller than the corresponding 
classical one TcRir) since (1 — f3)~^rc is very large. As 
a consequence, then, the corrected deflection angle turns 
out to be lower than the classical one. 

Inserting Ea. H15|l into Ea. p2|l . the lens equation may 
be conveniently written as : 



l-AA(/3,^,)(- 



P-2 



i?2 - ^s^ - - = 
2 



(16) 



^ Both 9 and 9s are measured in angular units and could be rede- 
fined 3.S 9 = £,/ Di and 9s = f^/Ds with ^ and tj in linear units. 

® For /3 outside this range, the deflection angle may still be com- 
puted, but there is not an analytical expression. 



"^ Hereafter, we will denote with the subscript GR quantities evalu- 
ated for n = 1, while we use the subscript FOG for the correction 
due to the fourth order theory. With this notation, for instance, 
the total deflection angle is a = csG-R + cepOG- Moreover, we 
refer to the results obtained for n = 1 as classical. 



with 



A^(/3,^c) 



V^(l - /3)r(l - /3/2) 
4^2r(3/2 - 13/2) 



(17) 



having defined d — O/0e.gr, with 0e,gb. the Einstein 
angle in the general relativity case given by : 



^E.GR — 



AGmDu 



(18) 



Note that, for n = 1 (/3 = 0), the corrected lens equation 
(fTB|l reduces to the classical one : 



d' 



,9,,9- 1 = 



thus ensuring the self consistency of the theory. As a 
general remark, the modified lens equation differs from 
the classical one because of the second term in the square 
brackets which furnishes a correction depending on both 
f3 and Vc (through 'dc = rc/Di9E,GR)- Actually, the mag- 
nitude of the correction strongly depends not only on the 
scaling radius r^ but also on the lens mass and the dis- 
tances {Di,Ds,Dis) since they all enter '&c- 



V. MICROLENSING OBSERVABLES 

The pointlike lens equation (|16|) differs from the stan- 
dard general relativistic one for the second term in the 
square parenthesis. Should this term be negligible, all 
the usual results of gravitational lensing are recovered. 
It is therefore interesting to investigate in detail how the 
corrective term affects the estimate of observable related 
quantities since, should they come out to be detectable, 
they could represent a signature of i?" gravity. 

Since we are considering the pointlike lens, a typical 
lensing system is represented by a compact object (both 
visible or not) in the Galaxy halo acting as a lens on 
the light rays coming from a stellar source in an external 
galaxy, like the Magellanic Clouds (LMC and SMC) or 
Andromeda. It is easy to show that, in such a configura- 
tion, the standard Einstein angle 9e,gr. and the images 
angular separation are of the order of few x 10^^ arcsec 
so that we are in the regime known as microlensing [26| . 

In the following, we will consider different observable 
quantities in microlensing applications in order to see 
whether it is possible to detect deviations from the stan- 
dard theory of gravitational lensing. As a fiducial system, 
we take a stellar lens in the galactic halo (with m = \ M© 
and Di = 20 kpc) and a source star in the LMC (with 
Ds = 50 kpc). The main results may be easily scaled to 
other values of (m, Di, Ds) by the corresponding Newto- 
nian Einstein angle. 

The dependence on the theory is parametrized by the 
two parameters (n, r^) entering the modified gravita- 
tional potential \i'A\ . Actually, we prefer to use (3 rather 



than n as first parameter since its value is limited in 
the narrow range (0, 1) with (3 = Q corresponding to the 
General Relativity case n — 1. In order to estimate a 
reasonable range for re, we proceed as follows. Let us 
rewrite the potential as : 



$ 



Gm 



-V 



1 



ri' 



P\ 



%- (1 



ri' 



P\ 



with rj = r/vc and $i = Gni/rc- Since <&! has the di- 
mension of a squared velocity, we introduce a further 
parameter vi 
reads 



1/2 

^^ SO that the potential conveniently 



$ = _- 



(1 



Comparing the two expressions for the potential, we get : 



Gm 



(19) 



with vi — vi/c. Rather than Vc, we will use in the follow- 
ing log-Di as parameter since its range is easier to eval- 
uate. Let us consider, for instance, the Sun as the field 
source and a planet as a test particle. Since the Earth 
velocity is ~ 30 km/s, we can indeed take (—6, —4) as a 
reasonable range for log-Di. As a final remark, note that 
Ea. (|19|) clearly shows why Tc depends on both the mass 
and the scale of the system under examination. 



A. The Einstein angle and the position of the 
images 

As a first application, let us consider the solution of the 
lens equation starting from the case of perfect alignment 
among observer, lens and source. Setting ds = 0, Ec ljlfcil) 
reduces to : 



l-AA(/3,^,)(- 



/3-2' 



^' = ; 



(20) 



which, for (3 = (n = 1) reduces to : 



i3' 



1 



^E.GR. 



i.e. the image is the well known Einstein ring with angu- 
lar radius equal to the Einstein angle defined in Ea. p8|) . 
In the general case, Ea. (|20(l must be solved numerically 
for given values of the model parameters (/3, log-Di). As 
a general result, the Einstein ring still forms, but the 
corrected Einstein angle 9e deviates from the classical 
one, i.e. -Oe 7^ 1. It is therefore interesting to study the 
magnitude of such deviations. To this aim, we define : 



£e 



7E,GR. 



df 



^E,GR 




FIG. 1: Contours of equal ee in the plane (/3, logwi) for 
£E = (-0.25,-0.20,-0.10,-0.05,0.05,0.10,0.20,0.25) from 
the bottom curve to the top one. 



which quantifies the relative deviation with respect to 
the classical result. Contours of equal Ee in the plane 
(/3,logwi) are shown in Fig. ^ 

Some interesting remarks may be drawn from this plot. 
As a general result, Ee could take both positive and neg- 
ative values, i.e. i?" gravity may increase or decrease the 
Einstein angle with respect to the classical result. Actu- 
ally, over almost the full parameter space, Ee is negative 
and, in particular, taking (/3,log'Di) ~ (0.58, —5) as fidu- 
cial values*, we find Ee — —0.25 which leaves open the 
possibility to detect the corrections due to i?" gravity. 

Before discussing this issue, it is worth examining how 
Ee depends on the model parameters. For fixed values 
of /3, the relative deviation is clearly an increasing func- 
tion of log vi , a result that can be easily explained qual- 
itatively. Indeed, the larger is logtii, the smaller is re 
and hence the earlier the correction term into the modi- 
fied potential starts contributing thus explaining why the 
Einstein angle deviates more and more from the classical 
one as log vi increase. A similar discussion may be done 
to explain why ee increases (in absolute value) with f3 
for fixed log vi . To this aim, let us remember that the 
deflection angle a changes from aoR for /3 = to aGR/2 
for /3 = 1 so that the lensing strength decreases with (3. 
As a consequence, the Einstein angle becomes smaller as 
P gets higher and hence Ee takes more negative values. 



Ind eed, fitting the LSB rotation curves, one gets (3 = 0.58 ± 0.15 
|20| . while log-Di ~ — 5 is suggested by considering a Jovian -like 
planet as test particle in the potential of a Sun -like lens. 



Modifying the Einstein angle has a deep impact on the 
microlcnsing phenomenology. As a first interesting appli- 
cation, let us remember that the duration of a microlcns- 
ing event is estimated as tE = Re/v±_ with Re = Di9e 
the Einstein radius and v±^ the transverse velocity of the 
lens with respect to the line of sight. Supposing that the 
distances {Di,Ds) and the transverse velocity v±^ were 
known (for instance, from deviations with respect to the 
standard lightcurve), a measurement of the Einstein time 
tE.obs translates into an estimate of the mass m of the 
lens. Denoting with /i the true mass (in units of 1 Mq) 
and with (J-gr the one estimated from using incorrectly 
the classical expression for the Einstein angle, it is re- 
spectively : 

tE,obs = fi^^'^tE{m = 1 Mq) = fJ-GRtE,GR{m = 1 Mq) 
so that we get : 



A* = 



t^GR 



d%{m 



1 Mr. 



(21) 



As can be inferred from Fig. ^ d\{m — 1 M0) < 1 over 
almost the full parameter space (/3,logui) so that the 
true lens mass is higher than the one estimated by the 
classical method. However, in most of the applications, 
one may statistically infer only the lens transverse veloc- 
ity, but not its distance so that what is indeed underesti- 
mated by the classical method is the quantity ^x{l — x) 
with X = Di/Ds- As such, therefore, taking into account 
the correction to the Einstein angle may lead to an in- 
crease of either the estimated lens mass or of the distance 
ratio X. Disentangling these two effects is quite difficult 
and need for non standard microlcnsing events. A pos- 
sible way out of this problem could be resorting to non 
standard microlcnsing events, but, in such a case, the du- 
ration of the event may differ from tE so that the above 
qualitative argument does not apply anymore. 

A similar discussion can be made for the optical depth 
T that measures the probability of finding a lens in the 
microlensing tube, a cylinder with main axis along the 
line of sight and with radius equal to the Einstein radius. 
Since the corrected Einstein radius is smaller than the 
classical one, the microlensing tube is narrower and hence 
the optical depth, for a given spatial distribution of the 
lenses, is lower than the classical one. 

As a second application, let us study how the images 
position change. In the Einsteinian case, the lens equa- 
tion may be solved analytically and one gets two images 
with positions given by : 



1? 



±,GR 



I {^s±yM 



+ 4 



(22) 



In the current case, the lens equation (|16|l must be solved 
numerically for given values of the source position -ds 
and the f{R) parameters {(3,\ogVi). We still get two 
images on opposite sides of the lens with one image lying 



inside and the other one outside the Einstein ring. The 
geometric configuration is therefore the same as in the 
standard case, but the positions are shghtly changed. To 
quantify this effect, we have studied the contours of equal 
e± in the plane (/?, log vi) with e± defined in a similar way 
as ee- Not surprisingly, we get the same results as for 
the Einstein angle, i.e. the contours are parallel to those 
of Ee but shifted by an amount depending on "dg- This 
is easily explained noting that the lens equation is quite 
similar to the one for the Einstein angle, while the term 
depending on the i?" gravity parameters is the same. 



B. Images amplification and Paczynski curve 

Gravitational lenses work also as a sort of natural tele- 
scopes amplifying the luminosity of the source in such a 
way to make visible objects that are otherwise too faint 
to be detected. As a general rule, the amplification is 
given by the inverse of the Jacobian of lens mapping. In 
the case of a lens with cylindrical symmetry (and hence 
also for the pointlike lens), this reduces to |23| : 



images. This result is not unexpected since the ampli- 
fication is defined in terms of the inverse of the Jaco- 
bian mapping and therefore increases when this latter 
decreases. Since det J is a decreasing function of /3, the 
behaviour of A is obvious. 

As a first important check, we have verified that the 
corrections to the amplification do not change the struc- 
ture of the critical curves defined as the loci in the lens 
plane where the amplification gets formally infinite. In- 
deed, using Eq. 1(201), we see that 6(^e) = 1 — a{'&E) that, 
inserted into Eg. 1)24(1 . gives det J (■(?£;) = and hence 
^{•dE) = oo. Therefore, the critical curve is still the 
Einstein angle as in the classical case and the only caus- 
tic (which is obtained projecting the critical curve on the 
source plane) is z?s = as usual. 

Since in a microlensing event the two images are 
merged together (given the tiny angular separation of 
the order twice the Einstein angle), what we observe is 
the total luminosity which is obtained by multiplying the 
source luminosity by the sum of the amplification of the 
two images. In the Einsteinian case, inserting Eas. H22|) 
into the right hand side of Eq. l(77|l , one gets for the total 
amplification : 



A 



'd d§ 



(23) 



Starting from Eq. 116(1 , after some algebra, one gets : 



A=\l-a^ -[{2- I3)a+{1- P)b + /3] b\ 



with : 



(2i92)- 
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(24) 



(25) 



(26) 



Note that we have considered the unsigned amplification 
since we are not interested to the parity of the images, 
but only to the luminosity variation. As a check of the 
consistency of the theory, it is easy to see that : 



\im A{^,P)=Agr = 



1-^ 



(27) 



so that the General Relativity result is recovered for n 
1. At the other extreme, we get : 



lim A(i}, (3) 



1 

4^ 



(28) 



so that the amplification is a increasing function of the 
slope parameter (3. This is in contrast with the results 
obtained for the Einstein angle and the position of the 



A^ 



tot,GR 



= A+ + A- = 
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(29) 



A similar relation for the corrected total amplification 
cannot be written since we do not have an analytical 
expression for the images position as a function of the 
source position. Actually, since the intrinsic source flux 
is unknown, we are unable to measure Atot so that detect- 
ing deviations from the standard theory is not possible 
at all. Nevertheless, a different strategy could be imple- 
mented. Assuming that the lens crosses the line of sight 
moving with a constant velocity, the source position as 
function of time t could be written as |23 : 
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(30) 



with -iJo = '9s{t — to), to the time of closest approach 
to the lens to the line of sight and tE the Einstein time 
defined above. As the time goes on, the source position 
changes according to Ea. (|Sn|) and hence the images po- 
sition varies. As a consequence, also the total amplifica- 
tion Atot becomes a function of time. In the Einsteinian 
case, inserting Eg. 1(30(1 into Eg. 1(29(1 . one gets the well 
known Paczynski lightcurve ,27j having the peculiar sig- 
natures of uniqueness and symmetry of the bump and 
achromaticity. Using the numerical solutions of Eg. ((16(1 . 
we can work out a corrected Paczynski lightcurve for the 
case of power -law f{R) theories. Some examples are 
shown in Fig.[21where, without loss of generality, we have 
set io = and tE = 10 d. 

As it is apparent from the figure, both the sign and the 
magnitude of the deviations from the standard lightcurve 
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FIG. 2: Paczynski lightcurve for Einsteinian gravitational 
lensing (solid line) and the R" gravity with (3 — 0.25 (short 
dashed) and (3 = 0.75 (long dashed). Upper (lower) pan- 
els refer to i9o ~ 0.1 (0.5), while left (right) panels are for 
logui = -5.5 (-4.5). 



FIG. 3: Relative deviation with respect to the Paczynski 
lightcurve of the corrected microlensing lightcurve for R" 
gravity with l3 = 0.43 (short dashed), /3 — 0.58 (solid) and 
/3 = 0.73 (long dashed). i9o ad log^i are set as in Fig. |3 



depend in a complicated way on both the theory pa- 
rameters (/Sjlogwi) and the minimum impact parame- 
ter ??o- Although a general rule cannot be extracted 
from the plots, it is interesting to note, however, that, 
for some combination of the parameters, the microlens- 
ing lightcurve may be also distorted with respect to the 
Paczynski one. Nevertheless, the typical signature of mi- 
crolensing events (symmetry, uniqueness of the bump and 
acromaticity) are preserved. 

To estimate quantitatively the deviations from the 
Paczynski lightcurve, we plot in Fig. |31 the quantity 
AAtot — 1 as function of the time t for the same val- 
ues of {t(),tE)- Considering, for instance, solid lines in 
the right panels, it is clear that the deviations are of or- 
der 10% which leaves open the possibility for a detection 
in the highest magnification events where the measure- 
ment errors on the lightcurve points are sufficiently low. 
A more detailed investigation is, however, needed to un- 
derstand whether such deviations may be mimicked by 
other astrophysical effects. Indeed, we have here con- 
sidered the case of a point mass acting as a lens on a 
point source. Should the lens (or the source) be a double 
system, this treatment breaks down and the curve may 
be distorted in a predictable way. Comparing these ef- 
fects with those introduced by _R" gravity is outside our 
aim here, but it is a task worth to be addressed in detail 
to understand whether there are detectable signature of 
power -law f{R) theories in the events lightcurves. 



VI. CONCLUSIONS 

Fourth order gravity has been recently proposed as vi- 
able alternatives to quintessence scalar fields and exotic 
fluids to solve the problem of cosmic speed up. In particu- 
lar, power- law f{R) theories have been shown to success- 
fully fit the SNela data without violating the constraints 



on the PPN parameters. Moreover, in the weak field 
limit, they give rise to a modified gravitational poten- 
tial that makes it possible to fit the LSB rotation curves 
without the need of any dark matter halo. 

Having successfully passed this impressive set of ob- 
servational tests, the proposed modification to Einstein 
general relativity worths to be further investigated con- 
sidering its effects on gravitational lensing. As a first 
step, we have derived an analytical expression for the 
deflection angle of a pointlike lens since this is the basic 
ingredient for a generalization to the case of extended sys- 
tems (such as galaxies and clusters of galaxies). Because 
of the deviations of the gravitational potential from the 
Newtonian one, the deflection angle turns out to differ 
from the Einsteinian one because of an additive power - 
law term depending on the slope P and the scalelenght Vc 
of the modified potential. For < (3 < 1, the corrected 
deflection angle is smaller than the classical one, but the 
amount of the corrections critically depends not only on 
/3 and Tc, but also on the lens mass m and the distances 
{Di,Ds,Dis). As a consequence of the corrected deflec- 
tion angle, the lens equation may no more be solved ana- 
lytically and the resort to numerical techniques is needed. 
The number of images is still two with an angular sepa- 
ration of the order of the corrected Einstein angle. 

To quantitatively investigate the impact of the correc- 
tions on observable quantities, we have considered the 
case of a solar mass object in the Galactic halo acting as 
a lens on a source star in the LMC so that we are in the 
typical conditions of microlensing regime. Denoting with 
£e the relative deviation of the corrected Einstein angle 
with respect to the classical one, we find out that both 
its sign and magnitude depend on the region in parame- 
ter plane (/3,logwi) considered, but significant deviations 
{ee ~ —20%) may be obtained for the value oi f3 (~ 0.58) 
suggested by the fit to the LSB galaxies rotation curves. 
A similar conclusion also holds for the deviations from 
the classical Paczynski lightcurve. Although a more de- 



tailed analysis (also including non standard microlensing 
events and observational uncertainties) is needed, these 
preliminary results suggest the intriguing possibility to 
detect i?" gravity signatures through a careful examina- 
tion of galactic microlensing. This may open the way 
to a completely new regime in which looking for correc- 
tions to the Einsteinian general relativity that may be 
complementary to the other cosmological probes. 

It is worth stressing that, since significant deviations 
from the standard theory of microlensing are possible, the 
current estimates of the optical depth r towards LMC (or 
any other target such as the galactic bulge, SMC or An- 
dromeda) and of the mean lens mass inferred from the 
time duration distribution of microlensing events could 
be altered. In particular, the classical estimates of the 
optical depth (the mean lens mass) should be revised 
downward (upward) for fiducial values of the parameters 
(/3,logwi). Interestingly, microlensing surveys towards 
LMC have concluded that less than 20% of the puta- 
tive dark halo is made out of compact baryonic objects 
(see, e.g., |29j and references therein) and there are mod- 
els explaining the observed optical depth as a result of 
known stellar populations only belonging to our Galaxy 
[33 and/or to the LMC [sj. On the other hand, our 
modified gravity makes it possible to fit the Milky Way 
rotation curve without any dark matter halo. Indeed, 
one may speculate that self-lensing models should be 
preferred in the context of power- law f{R) theories since 
one should have the possibility of explaining microlensing 
events without any dark matter halo (which is not needed 
in the f{R) theory we are considering). As a test, one 



could try to fit the Milky Way and LMC rotation curves 
using the modified gravitational potential and the mass 
distributions of the above quoted models. A successful 
result should be a serious evidence against dark matter 
halo in our Galaxy. 

The pointlike lens is only the first and essential step 
towards the formulation of a theory of gravitational lens- 
ing in a power -law f{R) theory. As a second step, we 
have to investigate the case of an extended lens such as 
a galaxy. This will offer us the possibility to get strong 
constraints on the theory parameters (/3, log f;i) by fitting 
to the images configuration in multiply imaged quasars. 
This will also represent a critical test of the f{R) theory 
since it is expected that /3 (and hence n) is the same for 
all the galaxy lenses considered. With more than 90 sys- 
tems (see, e.g., the CASTLES survey website |3^), this 
test will be statistically meaningful and will be presented 
in a forthcoming paper. 

As a final remark, we briefly comment upon the pos- 
sibility of evaluating the corrected deflection angle for 
other classes of f{R) theories. In principle, one should 
only insert the corresponding expression for the modified 
gravitational potential $ in the general equation ifTUI) for 
the deflection angle. Given the fourth order degree of the 
field equations, it is likely that the solution for $ may 
still be approximated (at least, within a limited range) 
by Ea. (|13|l with {f3, \ogvi) related to the f{R) parame- 
ters. It is therefore likely that the results discussed here 
for i?" gravity may be qualitatively extended to a more 
general class of higher order theories of gravity. 
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